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Q ■ ABSTRACT 

. We propose two methods that enable us to obtain approximate solutions of the lens 

equation near a fold caustic with an arbitrary degree of accuracy. We obtain "post- 
linear" corrections to the well-known formula in the linear caustic approximation for 
the total amplification of two critical images of a point source. In this case, in order 
to obtain the nontrivial corrections we had to go beyond the approximation orders 
' earlier used by Keeton et al. and to take into account the Taylor expansion of the lens 

, equation near caustic up to the fourth order. Corresponding analytical expressions 

are derived for the amplification in cases of the Gaussian and power-law extended 
^ ^, source models; the amplifications depend on three additional fitting parameters. Con- 

ditions of neglecting the correction terms are analysed. The modified formula for the 
amplification is applied to the fitting of light curves of the Q2237-I-0305 gravitational 
^ ' lens system in a vicinity of the high amplification events (HAEs). We show that the 

, introduction of some "post-linear" corrections reduces by 30 per cent in the case of 

known HAE on the light curve of image C (1999). These corrections can be important 
, for a precise comparison of different source models with regard for observational data. 

■ Key words: gravitational lensing: micro ~ quasars: individual (Q2237-I-0305) - grav- 

I itational lensing: strong - methods: analytical 

a^ ' 



1 INTRODUCTION 

An extragalactic gravitational lens system (GLS) forms several images of a single quasar. The light from the quasar intersects 
• • , a lensing galaxy in different regions which correspond to different images. Variations of gravitational fields in these regions due 
_ ^ ' to stellar motions are practically independent and lead to independent brightness variations in different images (gravitational 
I microlensing) . Comparison of the l ight curves of different images a llows one to obtain a valuable information about the lens 
$H ' itself and about the source as well ( Schneider. Ehlers fc Falkolll992l : IWambsgansl I l200d) . One of the important applications of 
5^ , this effect deals with a unique possibility to study a fine structure of the central quasar region with the use of GLS. This idea 



tnis ettect deals witn a unique possibility to study a nne structure ot tne central quasar region witn tne use ot (j,Ls>. inis idea 
first proposed by Grieeer. Kavser 8z Refsdal (j 19881 ) appeals to the high amplification events (HAEs) in one of the images o f 



the quasar in GLS. An interesting applications of HAEs are known also in Galact ic microlensing (see, e.g. Wambseans j 2006 ) 



2004 ) 



in pa rticular, using the caustic crossing events for resolving of stellar profiles IjBogdanov fc Cherepashchukll2002l : I DominikI 



The conventional explanation of HAE relates it to the caustic fiel d in the source plane f ormed due to the inhomogeneous 
gravitational field of a lensing galaxy on the line of sight of the image ([Schneider et al.ll 19921 ). The source crossing of a caustic 



leads to a considerable enhancement of the image brightness, the crossing of a fold caustic being the most probable. The 
corresponding variations of the brightness in a neighborhood of HAE can be approximately described by a formula containing 
a fe w fitting parameters . This makes it possible to estimate certain HAE characteristics, in particular, such as the source 



size ( Grieger et al. 19881 ) . For example, in the case of the well-known Q2237-I-0305 GLS (Einstein Cross), several HAEs was 
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observed ( Wozniak et al Iboool : [Alcalde et alJboodJudalski et alJbooel). and the estimates of the source size have been obta ned 
within different source models Iwvithe. Webster fc Turneilll999l . I2OO0I : IWvithe et allboool : lYoneharallioOll : IShalvapinllioOl 



Shalvapin et alJbood : iBogdanov fc Cherepashchukl2002). Almost all HA Es in the Q2237+0305 GLS are attributed just to a 



fold caustic crossing i n the source plane (see, e.g., ICil- Merino et al.ll2006l l. A possibility to distinguish different source models 
is also discussed ( e.g. Goicoechea et al. 20031 ') . 

The lens equation near a fold can be expanded in powers of local coordinates; in the lowest orders of this expansion, 
the caustic is represented by a straight line; so, this approximation is often referred as a "linear caustic approximation". In 
this approximation, the point source flu x amplification is given by a simple formula, which depends on the distance to the 
caustic and contains two parameters (e.g. Schneider et al. 19921 . Cassan 2008). In most cases, the linear caustic approximation 
is sufficient to fit the observed light curves over the range of HAEs at the modern accuracy of fiux measurements. The need 
for a modific ation of this formula - e.g., by taking the caustic curvature into ac count - is, nevertheless, being discussed for 
a long time jPluke fc Websteij|l999l : IShalvapinI [20011 : IPeicha fc Hevrovskvl [2OO9I ') . We hope that the future improvement of 
the photometric accuracy will make it possible to obtain additional parameters of the lens mapping, which are connected 
with the mass distribution in the lensing galaxy. At the same time, we will show that the consideration of "post-linear" 
terms is sometimes appropriate to explain the present observational data. We note that the corrections to the amplification 
factor in the case of macrolens ing were the subject of investigations dealing with the problem of "anomalous fiux ratios" 

i Keeton. Gaudi fc Petteilboosi) . ,— — , , 

Since the work WlKochanek' ("2004') followed b y a number of works llvlortonson, Schechter fc W ambsgansj 2005': 'Gil-Merino et al. 



s mce tne worK p yiivocnaneif ( zuu4 ) ii 

20061 : IVakulik et ahl EoOT: Anguita ct al. 



2008; 



Poi ndexter. Morgan fc Koc hanck 2008; Poindcxtcr fc Kocha nck 2010a 
merous statistical methods have been developed to process complete light curves. This approach is very attractive because 
it allows one to consider the whole aggregate of observational data on image brightness variations in order to estimate the 
microlens masses and the source model parameters. However, this treatment involves a large number of realizations of the 
microlensing field, which requires a considerable computer time. On the other hand, the source structure manifests itself only 
in HAEs; far from the caustic, the source looks like a point one, and all the information about its structure is lost. If we 
restrict ourselves to the HAE neighborhood, then we use the most general model concerning a microlensing field described by 
a small collection of Taylor coefficients in the lens mapping. Therefore, the low-time-consuming semi-analytical investigations 
dealing with caustic neighborhoods still preserve their importance, not speaking about their use in computer codes. 

In this paper, we propose relations for the total magnification of two critical images of a point source in the first nontrivial 
"post-linear" approximations and use them to modify the magnification of an extended source near the fold caustic. In order 
to obtain non-zero corrections to the total amplification of two critical images, we have to consider additiona l higher-order 
terms in the expansion of the lens mapping in comparison with the earlier works (see, e.g.. lKeeton et al.ll2005r ). Though the 
corrections are expected to be small, they appear to be noticeable in some cases even in an analysis of the existing data on 
light curves in the Q2237-I-0305 GLS. The structure of this paper is as follows: after the derivation of approximate solutions 
of the lens equation in the required order, we deduce a formula for the magnification of a point source. The result is used 
to obtain the magnification of a small Gaussian source near the fold caustic. We use the formula for the amplification of 
an extended source to fit the light curves in GLS Q2237-I-0305. The obtained post-linear corrections appear to improve the 
agreement with observational data. 

Appendices |A] and |B] contain details on the approximation methods involved. Power-law source models are cons i dered in 
Appendix [C] Explicit comparison of our approximate formulas with counterparts from the paper bv lKeeton et al.l ( 20051 ) is 
made in Appendix [D] 



2 INITIAL EQUATIONS 



First, we recall some general notions concerning the gravitational lensing that can be found, e.g., in the book bv lSchneider et al 
( 1992h . The normalized lens equation has the form: 



y = x-V^{x), (1) 

where $ {x) is the lens potential. This equation relates every point x = (xi, X2) of the image plane to the point y = (j/i, 2/2) 
of the source plane. In the general case, there are several solutions X^i) (y) of the lens equation ^ that represent images of 
one point source at y; we denote the solution number by the index in parentheses. 

If there is no continuous matter on the line of sight, the potential must be a harmonic function A<1> — 0. Below, we will 
assume that this condition is fulfilled in a neighborhood of the critical point. We note, however, that if the continuous matter 
density is supposedly constant during HAE, this can be taken into account by a suitable rescaling of the variables. 

The amplification of a separate image of a point source is 

7^(,)(y) = l/|j(X(,)(y))|, (2) 
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where J (x) = \D{y)/D{x)\ is the Jacobian of the lens mapping. In the microlensing processes, microimages cannot be 
observed separately; therefore, we need the total amplification which is a sum of the amplification coefficients of all the 
images. 

The critical curves of the lens mapping ([T|) are determined by the equation J (cc) — and are mapped onto the caustics 
in the source plane. The stable critical points of a two-dimensional mapping can be folds and cusps only, the folds being more 
probable in HAE. In this paper, we confine ourselves to the consideration of fold caustics. When a point source approaches 
the fold caustic from its convex side, two of its images approach the critical curve, and their magnification tends to infinity. 
They disappear when the source crosses the caustic. These two images are called critical. 

The standard consideration of the caustic crossing events deals with the Taylor expansion of the potential near some 
point Per of the critical curve in the image plane. Let this point be the coordinate origin. We suppose that ([T]) maps p^r onto 
the coordinate origin of the source plane. Further, we rotate synchronously the coordinate systems until the abscissa axis on 
the source plane be tangent to the caustic at the origin. The quantity |j/2| defines locally the distance to the caustic, and yi 
defines a displacement along the tangent. For the harmonic potential, we can write 

yi = 2xi + a (^Xi — a;!) + 2bxiX2 + c (^Xi — Sxix?,"^ — d — 3a;2xf) + gx^ + 

y2= h [xi — xl) — 2axiX2 + d (s? — Sxixfj +c(^xl — 3x2x1) + fx^ + (3) 

where a, 6, c, d, g, and / are expansion coefficients. If the y2 axis is directed toward the convexity of the caustic, then 6 < 
(at fold points, 6 7^ 0). 



3 EXPANSION OF THE CRITICAL SOLUTIONS IN POWERS OF A SMALL PARAMETER 
3.1 Method 1 

We now proceed to the derivation of approximate solutions of Eqs. Q. To do this, we present two different methods which will 
be used to have a possibility of mutual checks of cumbersome calculations. The first method deals with analytical expansions 
in powers of a small parameter. However, it results in nonanalytic functions of coordinates leading to nonintegrable terms in 
the amplification factor. The second method does not lead to such problems, though it uses a somewhat more complicated 
representation of the solution of the lens equation (containing square roots of analytic functions). The methods agree with each 
other in a common domain of validity; moreover, we use the second method to justify some expressions in the amplification 
formulas in terms of distributions to validate appl ications to extended source models . 

First, we use a regular procedure proposed by Alexandrov. Zhdanov &: Fedoroval ( 2003h to construct solutions of Eqs. ^ 



with a desired accuracy. This procedure is useful to study the light curve of the point source which has a trajectory crossing 
the fold caustic at some nonzero angle. We suppose that the source and the caustic lie on different sides from the y\ axis. 
Then, for j/2 > 0, we substitute 

yi = t^iji, xi = t^xi, X2 = tX2, (4) 

where i = 1,2, and t can be considered as a parameter of vicinity to the caustic. This is a formal substitution that makes 
easier operations with different orders of the expansion. After performing calculations, we shall put t = 1 and thus return to 
the initial variables yi. However, if we put yi to be constant with varying t, then this substitution allows us to study a local 
behaviour of critical image trajectories; t = corresponds to crossing the caustic by a point sou rce, and can be consi dered 



as the time counted from the moment, when two critical images appear. Indeed, one can show ( Alexandrov et al. 20031 ) that 



two critical solutions Xi{t) = (X(i)i, X(i)2) of the lens equation can be represented by analytic functions of t that, in the above 
special coordinate system, have the behaviour Xf^ji oc t^, X(i)2 oc t (see Appendixl^. This allows us to look for solutions of 
Eqs. ((3|, by using the expansions of Xi in powers of t: 

Xl = Xio + Xlit + X12 + 

X2 = X20 + X2lt + X22t^ + ■■■■ (5) 

It should be stressed that the analyticity in t does not mean that the coefficients of expansions ([5]) will be analytic functions 
of coordinates iji in the source plane (see below). 

In terms of the new variables (|4]), system Q takes the form (up to the terms ~ t^) 

y\ = 2x\ — ax2 + t{ 2hx\X2 — dip) + {ax\ — 3cxiX2 + 9X2) , 

y2 = —bx2 + t{~2ax\X2 + cx'2) + t^{bxi — 3dxiX2 + fx2). (6) 

The substitution of expansions ((Sj into ([6| allows us to determine all coefficients successively. The results of calculations 
are as follows. 

For the zero-order terms: 
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X20 = £\/v2l\h\ 



where e = ±1 determines two difTerent solutions. 
The first-order terms are 



262 



^21 = - 



62 



-y2 - -pji 



(7) 



(8) 



where f?^ = + 



}? . The solutions up to this acc u racy h as been obtained earlier by Alexandrov et al. 1 2003[ ) and Keeton et al 



(see also ICongdon. Keeton fc Norderren I l|2008h ). The contributions of this order are cancelled in calculations of the 
total amplification factor of two critical images. Therefore, to obtain a nontrivial correction to the zero-order amplification, 
higher order approximations should be involved. 

The second-order terms contain an expression nonanalytical in j/2: 



2:12 



^ (3a^ + 5a^b^ + 206" - 26^d - 2b'^g + 3ac^ - 6a^c + 36cd - 8a^bd + 2abf) yl + 



263 



{2a^c - b^c + 3abd - 2a'^R^- b^R^) yiya + ^y?, 



1/2 

6 



(lOa^c - 5c^- Sa^i?^ + lOabd - 46/) §2 (ac + bd - aR^) yi - 

863 ^ ; » 452 V ' " 9,b y2 



(9) 



(10) 



3.2 Method 2 



The second approach to the construction of approximate solutions of the lens equation in a vicinity of the fold is described in 
Appendix [B] This allows us to provide the critical solutions of system ((6| in the following form: 



I + treJw, 



X2 



tS + £y/w, 



±1. 



(11) 



The iterative procedure described in Appendix iBl yields analytical expansions for the functions p,r,s,w both in powers 
of t and iji. The application of this method to system (|6]) gives (up to the terms ~ t^) 



Xio + t X12, 

where xio and X12 are given by relations (O and ((9]), and 



26 
6 



■ + bd)] 



s — yi 

2b" 



462 



6 



262 
1 



-J/2, 



R^yi + - (bd — ab^ + ac — a'^) yiy2 + {5a^(i?^ 



2c) + 5c^ - 10a6d + 4b f) 



J/2 



(12) 



(13) 



(14) 



As is seen, all these expressions are analytic functions of both t and yi. If we expand ^/w in powers of t, then we 
immediately have the solution in the form ((SJ with coefficients (O, ((S)), (©, and (fTO)) . 



4 AMPLIFICATION FACTOR 
4.1 Point source 

For the Jacobians J (-X'(i) (^'^y)) (* = 1, 2 corresponds to e = ±1), we obtain up to the terms ~ 



J = 4:ety/\b\ y2 + 4t' 



(R^-c)^ 



-y-i 



2/2 



62 



-J/1 



+ 7^ (7a^i?^ - 8c_R^ + Ic - 6a^c - 30a6d + 246^c + 126/) 
263 \ I 



J/2 



26 ^2 



(15) 



According t o (l2l). the value of yields the amplification of individual images. Note that the terms up to order ~ 1? 
were obtained by Keeton et al. ( 2005^ . The final result for the total amplification of two critical images (in terms of the initial 
variables yi after putting t = 1) is as follows: 



K 



1 e(y2) 



1 + Py2 + Qyi 



4 y2 



(16) 



where 
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P^2n+^[aW + {a^ -cf]- ^{2f -5ad), Q ^ ^ (^a' - ac + ab^ - bd) , k = , (17) 



is the Heaviside step function. Note that k is the caustic curvature at the origin jCaudi fc Petterdbood : lAlexandrov et al 



2003) which enters exphcitly into the amphfication formula. 



Formula (|16[) yields an effective approximation for the point source magnification near the coordinate origin provided 
that y2 > 0, and 1/2/2/1 is not too smaU (see the term containing a). For a fixed source position, this can be satisfied always 
by an appropriate choice of the coordinate origin, so that the source will be situated almost on a normal to the tangent to 
the caustic. 

If the source is on the caustic tangent or in the region between the caustic and the tangent, then formula (|16|) does not 
represent a good approximation to the point source magnification. Nevertheless, in case of an extended source, we will show 
that result (|16|) can be used to obtain approximations to the amplification of this source even as it intersects the caustic. 
However, to do this, we need to redefine correctly the convolution of (|16p with a brightness distribution. 

4.2 Transition to extended source 

Let I{y) be a surface brightness distribution of an extended source. If the source center is located at the point Y = {Yi, Y2) 
in the source plane, then the total microlensed flux from the source is 



F{Y) = JJ I{y{x) - Y) dx^dx2 = jj K{y)I{y - Y) dy,dy2, (18) 
where the point source amplification K{y) = "^Ki is the sum of amplifications of all the images. The result of using the 



first integral from Eq. (|18|l obviously is equivalent to the result of the well-known ray-tracing method ([Schneider et al.lll992l 'l 
(when the pixel sizes tend to zero). Near a caustic, one can approximate K{y) = Ko + Kcr{y), where Ko is an amplification 
of all noncritical images that is supposed to be constant during HAE, and Kc- is the amplification of the critical images. 

Formula (I16|) contains the non-integrable term ~ Q{y2){y2)~^^^ ■ Therefore, the question arises of how formula (|16[) can 
be used in situation when the extended source intersects a caustic and some part of the source is in the zone between the 



tangent and the caustic. In view of Section [3.21 it is evident that the mentioned term is a result of the expansion of the root 
\/y2 + K,yft^/2 + ... in the approximate solution (|11)I14|I . Any non-integrable terms in Kcr does not arise without using this 
expansion. It is easy to show that, in order to define Kcr correctly, one must replace the term 6(1/2) (2/2)"'^^'^ in (|16|) by the 



distribution (generalized function) (1/2)+^'''^ (iGel'fand fc Shilovlll964 ). We recall that the distribution y^'''^ of the variable y 
is defined by the expression 

00 00 



for any test function f{y). 

After this redefinition, we have 

Kcr = [1 + Py2 + Qyi]~ ^j=yl{y2)X'^\ (19) 

This formula can be used to correctly derive an approximate magnification of a sufficiently smooth extended source including 
the case where the source crosses the caustic. 

4.3 Extended Gaussian source 

Now we use formula H19p to derive the magnification of a Gaussian source with the brightness distribution 

/G(2/) = 4^exp(-yVi'), (20) 

where the parameter L characterizes the source size. 

The amplification of an extended source is defined as the ratio of the lensed fiux (|18p to the flux of the unlensed source 
Fo = JJ I (y) dy\dy2 which is equal to 1 in case of formula (|20|) . The amplification of a Gaussian source Kq is obtained by 
the substitution of ((20ll and gH) into (fTS)) . 

Further, we introduce the dimensionless coordinates s = Y\ / L, h — Y2/ L oi the source centre and the functions 

7 -1 °° r (i -I- '=+" '1 

h{h)= / n''-^/"exp(-it" + 2tt/t)rftt^ ^ V , ' ' {2hf. (21) 

/ ^ '2 ^ — ' n! 

11=0 



-3/2 
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Figure 1. Behaviour of individual functions 1)2411 which generate dependence ||23| I and the ratios 1)2611 of correction functions to the main 
one. 



These functions can be e xpressed in terms of the confluent hypergeometric function \F\ or the parabohc cyhnder function 
D ( Bateman fc Erdelvillliii '): 



h (h) = ir (1 + I) (1 + 1,1; h') + /^^ (I + I) iF, (| + |, |; /,^) = 2-(4 + i)r (fc + i) e^i3_(,^ . ) (-^/2 . h) 



The substitution of ^ and JgOj in (HHJ yields 



/^G (s, ft) 



Here, 



|$o (ft) 



P$l (ft) - -$2 (ft) + QS^O (ft) 



Ks2$2 (ft) 



} 



(22) 



(23) 



(24) 



$0 (ft) = lo (ft) exp (-ft") , $1 (ft) = Ji (ft) exp (-ft^) , $2 (ft) = [hlo (ft) - h (ft)] exp (-ft^) . 

We have checked this result by the direct substitution of expansions ((6| in the first integral of Eq. p8|) . assuming t — \fL 
and expanding the resulting integral in powers of this parameter up to the second order. Note that the main term of ()23|) 
which corresponds to the linear caustic approximation was first obtained by Schneider fc Wei3 (1987). 

Analogous considerations allowed us to obtain formulas for the magnification of extended sources for two types of power- 
law brightness profiles (Appendix [C]); the results are represented analytically in terms of hypergeometric functions. 

Let us discuss formula ()23D in more details. The functions $i (ft) are shown in Fig. [T] The distinctive variations of these 
functions take place for —21/ < j/2 < 2L. For j/2 < — SL, they are practically equal to zero, and, for j/2 > 2L, they have the 
following asymptotic behavior: 



(25) 



*o(ft) = ./T, $i(ft) = V7rft, 
V ft 

We also introduce the functions 
$3 (ft) = -J-i/^o, $4 (ft) = -*2/*o, (26) 

which are designed to estimate the contribution of the correction terms. These functions are also shown in Fig. [1] At the 
origin, we have $3(0) = $4(0) = 0.338. 

A special attention must be given to the term which is proportional to the monotonically increasing function $1 (ft). 
This correction becomes especially noticeable on the inner side of the caustic, its sign being determined by the sign of the 
parameter P from Eqs. ()17p . Such behaviour allows us to hope that the determination of this term from the observational 
data will not be too difficult. 

The effect of the third term of formula (|23|) that is proportional to "I>2 (ft) is shown in Fig. 2 by the dependence 
$0 (ft) — Q'1'2 (ft) for various values of the coefficient a. According to ((23}, the parameter a — Lk/2 equals to half of 
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Figure 2. Influence of the correction term proportional to <3>2- 



the ratio of the source radius to the curvature radius of the caustic at the origin. One can see from Fig. [2] that this term can 
noticeably affect the determination of the source size and the time moment of crossing the caustic . 

The last two terms of Eq (|23|) contain s cx j/i and cx j/i^. For a fixed position of the source, one can exclude these 
terms by the origin displacement. In case of the rectilinear motion of a source, the impact of these terms can be noticeable 
for small angles of the intersection of the source trajectory with the caustic. 

We now discuss the conditions of applicability of the approximation methods involved for modelling the light curves in 
a vicinity of HAE. For example, with regard for observational data, we can require that, in the interval —0.5 ^ /i ^ 1, the 
contribution of each correction do not exceed five per cent. We see from Fig. [T]that, in this interval, $3 (h) < $3 (1) ~ 0.83 
and |$4 {h)\ < $4 (—0.5) ~ 0.73. From the condition of smallness of the corrections, we find the following restrictions on the 
parameters of the model: 

LP < 0.06, Lk<0.13, LQ • cot (/3) < 0.05, Lk • cot^ < 0.27, (27) 

where /3 is the angle between the source trajectory and the caustic. Under these conditions and within the specified margin 
of error, one can use the linear caustic approximation. Thus, in addition to the requirement of the smallness of the caustic 
curvature, we obtain some extra restrictions on the possible variations of the lens potential on the scale of the source size. It 
is clear that the strengthening of requirements for the model accuracy and/or the extension of the interval of h leads to the 
strengthening of the conditions found. 

Clearly, formulas (I19|l and (|23p make sense only if we can ignore the discarded high-order terms. This means, in turn, 
that the correction terms in (|23|) themselves must be sufficiently small. 



5 HAE IN THE EINSTEIN CROSS 

The Einstein Cross QSO 2237+0305 I Huchra et al. 20081 ) consists of a quadruply imaged quasar and a lensing galaxy that 
is the nearest of all known gravitational lens systems {zg = 0.0395). The gravitational delay times between images in this 
GLS are of the order of h ours. This follows from a highly symmetric confi g uration of the images and is partially con firmed 



Fedorova et al. 20081'). Since 



(jfLb are ot tne order ot h ours, inis tollows trom a nignly symmetric conn g uration ot tne imag ( 
by observations (see, e.g., Schmidt. Webster fc Lewi3 19981 : Dai et al. 20031 : Vakulik et al. 20061 : 
the Einstein Cross is a very suitable object for microlensing studies, its images have been continuously monitored by different 
groups for more than a dozen of years. In this system, significant microlensing-induced brightness pe aks on light curves of the 
quasar images were detected (see, e.g., Woziiiak et al. 200ol : Alcalde et al. 2OO2I : Moreau et al. 20051 ). 
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Fitting the light curves and estimations of HAE parameters in GLS Q223T+0305 



We now apply formula (|23p to the fitting of the light curves near HAE. For a moving source, Yi = Vi{t — tc), where t is the 
time, tc is the time of the crossing of the caustic by the source centre, Vi is the projection of the source velocity on the axis yi. 
We suppose that V2 » Vi, i.e. the source crosses the caustic effectively and does not move along it. Our numerical simulations 
have shown that the terms depending upon the coordinate s contribute only for small angles between the source trajectory 
and the tangent to the caustic. Therefore, we do not take them into account, and, correspondingly, the parameter Q is not 
involved into consideration. Introducing the parameter T = L/| V2I, we obtain h — ±{t — tc)/T (the sign "+" corresponds to 
the source motion along the positive direction of the 1/2 axis). 

We consider the know n HAE in the light curve of image C of GLS Q2237+0305 using the OGLE data recorded during 
1999 jWozhiak et al.ll2000l ). Let Fq be the flux from image C when the microlensing is absent. Under the supposition that 



the proper brightness variations of the quasar in GLS can be neglected and taking expression H23[) for the amplification into 
account, we obtain the formula for fitting the fiux from the microlensed Gaussian source, 

f"' {t)^A + B$o (h) + C$1 (h) + D$2 (h) , (28) 
which contains the parameters 



A = FoKo, B = Fo/ y/inT\b\\V2\, C = FoPy/T \V2\/4n \b\, D = -^FoKy/T \V2\/n \b\ , 

and tc and T which appear nonlinearly. The quantity Ko in the expression for yl is a part of the amplification due to 
noncritical images. The parameters A and B are evidently positive, D is negative, and C can have values of both signs. 
As discussed above, the possibility to use the linear caustic approximation or formula (|23p is determined by the ratios of 
corrections coefficients to the coefficient B of the zeroth approximation: C/B = LP, D/B = —Lk/2. 
To fit the light curve, we used the minimization of the weighted sum of squares: 

N 

S = Y,W,[F-F^' {t,)]\ (29) 
1=1 

where Fi is the result of the i-th measurement, and Wi = l/crf is its weight that is expressed through the corresponding 



dispersion estimate ai (jWozhiak et al. 200Q '). The fitting quality is often characterized by the parameter x = Suii-n/v, v being 



the number of degrees of freedom. The value of this parameter in the optimal case should tend to 1. 
As F^'^ , we have considered the following models: 

f" (t)^A^ B^o (h) , 

F^ {t)^A + B^o (ft) + C$1 (h), 

F^ {t)=A + B$o (ft) + D'i'2 (ft) . 

We also analysed the model that takes both correction terms C$1 (ft) + D^2 (ft) into account. However, we found that it does 
not allow us to obtain the coefficients that are statistically significant simultaneously. For comparison, we also considered the 
model with a correction term line ar in ft; such a correction can be caused by the own variability of the quasar, or by the 
influence of noncritical images (cf. Yoneharal 2001 ) : 



F^ {t)=A + B^o (ft) + Eh. 

The results of best-fitting with different models are presented in Table [T] It contains the estimates of model parameters 
and their central 95-per-cent confidence intervals that have been found by the Monte-Carlo simulations under supposition of 
the normal distribution of errors. In all the models, the correction terms are statistically significant. The probability that the 
correction coefficient is occasionally nonzero is certainly less than 10"'^ in every case. On the other hand, all three models can 
compete with one another on an equal footing (and probably with another effects such as those due to a complicated source 
structure). 

Note that the fiux variation of the model light curve is roughly equal to 0.18 mjy, and the standard deviation of data is 
a ~ 0.006 mJy. Thus, the assumption of the 5-per-cent tolerance in a vicinity of the light curve maximum, which has been 
used in criteria (I27|l . is rather realistic. In the case of the F^-model, we find C/B = LP « 0.076. The comparison with the 
first inequality in (|27|) indicates the agreement with the statistical significance of the first correction term. 

Fig. [3] shows the results of light curve fitting for image C in the region of HAE with the _F^-model. Here, 44 data points 
were included corresponding to the [1289,1442] epoch interval (light circles). To choose a fitting interval in a vicinity of the 
light curve maximum, we used the fact that the values of fitting parameters are reasonably stable with respect to a reduction 
of the interval. This allowed us to find a preliminary estimate of the "background level" A. After that, all the points above A 
have been involved in the final treatment. 

The models F^ and F"^ fit the data equally well. Let us discuss the model F^ and the role of the second correction. As 
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Table 1. Model parameters an d their central 0.95 c onfidence intervals as the results of light curve fitting for Q2237+0305C. 
The data due to OGLE group l|Wozniak et al.ll200d 'l. 



Fitting tc 


T 


A 


B 


Correction Smin ^ 




model (JD-2450000) 


(days) 


(mJy) 


(mJy) 


coefficient (mJy) 





^2 

F3 



1384.7ti;o 



1382.6' 



-1.5 



1367. if |j 
1383.3tJ;4 



33.0 
34.7 
40.1 
34.9 



+2.7 
2.6 
+2.9 
2.7 
+3.2 
3.5 
+3.2 
2.6 



^+0.006 
-0.006 
"+0.006 
-0.006 
3+0.006 
-0.007 
n qt;a+0 007 

U.OOO_Q ggg 



0.364: 
0.366j 
0.363^ 



r. + 0.003 

-0.003 
3+0.004 
-0.004 
7+0.005 
' -0.004 



0.076: 
0.079^ 
0.067"^ 



C-- 
D : 
E -. 



-0.006^ 



f0.003 
-0.004 
-0 071 +0 01'' 

-0 004+° 0°2 



53.91 


40 


1.35 


38.43 


39 


0.99 


34.97 


39 


0.90 


39.38 


39 


1.01 



Notation for the fitting models and the parameters are given in the text. 




Figure 3. Observed light curve of Q2237-I-0305C (OGLE data) and the best-fitting with the model (solid line). 



compared with the other models, the model leads to a somewhat lower value of x^, but the confidence interval for tc is 
four times wider than that in the case of the model F", and the source size appears to be 20 per cent more than that within 
other models. Here, the typical features of the second correction mentioned in the previous section become apparent. The true 
determination of D from formula (|28p could allow us to estimate the ratio of the curvature and source radii: nL = —2D/B. 
But the obtained value of D corresponds to the curvature radius Rc = l/n that is about twice less than the source "size" L. 
So, the condition of smallness of the correction is violated, and the model F^ is not acceptable. On the other hand, neglecting 
the second correction in (|28[l (i.e., the transition to the model F^) corresponds to Rc » L. Our calculations have shown 
that the a priori assignment Rc > L and the introduction of the appropriate additional term into the model F^ have little 
influence on the estimates of A, B, and C. 

To sum up, the model F^ fits satisfactorily the data on HAE in question within the present accuracy, though we cannot 
rule out competing models and/or competing effects. 

We also analysed HAE in the light curve of image Q2237-I-0305A that happened in 1999 jWoziuak et alibood ). However, 
the treatment of this HAE does not allowed us to find statistically significant estimates of the corrections because of lacking 
the data corresponding to the caustic inner region. 



6 DISCUSSION 

We proposed two methods that enabled us to obtain the critical solutions of the gravitational lens equation near a fold caustic 
with any desired accuracy. The first method is based on the expansion of the solutions and the amplification factor Kcr of 
two critical images of a point source in powers of a formal parameter t that can be interpreted as a parameter describing the 
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proximity to the caustic. We derived a representation of Kcr which is convenient to obtain the amplification factor for an 
extended source. This representation is not analytic in local coordinates, and Kcr contains a nonintegrable term. In order to 
elucidate this term, we considered a different representation of the solutions which contains analytic functions of coordinates 
and the parameter t, as well as square roots of such functions. This method yields only integrable expressions and explains 
the meaning of the nonintegrable term in Kcr'- it should be considered as a distribution which can be used to calculate the 
amplification of a small (and smooth) extended source near the fold caustic by taking the appropriate convolution. 

In order to obtai n nontrivial corrections t o Kc r, we had to use hig her orders of the expansion of the lens equation as 
compared to works bv I Alexandrov et al. I (|2003l ) and lKeeton et al.l (|2005l ). This is a consequence of the cancellation of terms 
~ O (t^) which are present in the amplifications of separate critical images. Therefore, we had to deal with terms of the order 
~ O (t^^ and the Taylor expansion including the 4-th order in the lens equation ([3]). The modified formula for Kcr contains 
3 extra parameters that are combinations of 5 coefficients of the Taylor expansion (we confined ourselves to the case of no 
continuous matter on the line of sight). 

Based on our result concerning Kcr, we derived an asymptotic formula for the amplification Kq of a small Gaussian 
source. We analysed the role of new "post-linear" corrections and formulated the conditions of applicability of the linear 
caustic approximation. The fitting of the light curve of GLS Q2237+0305C shows that some of these corrections can be 
statistically significant. 

Of course, we are far from the thought that the available observational data make it possible to determine all additional 
parameters simultaneously (and even more so, the coefficients of the expansion of the lens mapping up to the orders involved). 
Such a determination requires a considerable improvement of the photometric accuracy in the future. Our result is more 
modest: we have shown that, even at the present level of accuracy, some elements of the light curves require to consider the 
higher-order terms in the solutions of the lens equation. At least, after the introduction of some post-linear corrections, 
decreases by 30 per cent thus approaching 1. One may hope that an enhancement of the observational accuracy will increase 
the role of post-linear approximations. 

In this connectio n, we note that, b esides estimat ing the so urce size in GLS on the ba sis of light curves, a number of 
authors Ishalva pin 2001 Shalvapin et al. 200 2: Bo gdanov fc Cherepashchuk 2002 : Goicoc chea et al. 20031 : iKoch anck 2 0o3 : 



Mortonson et al.1 [2005 ; Gil-Merino et al. 20061 : Ivakulik et al, 2007 ; Anguita'et al l feoo^ ) discussed som e delicate question s 
concerning a fine quasar structure. For example, Goicoechea et aL I 2003h wro te that t he GLITP data ( Alcalde et al.l [20021 ) 
on GLS Q2237-I-0305 admit only accretion disc models fsee also iGil-Merino et al. 2006: lAnguita et al. 2008l l. Obviouslv. the 
presence of an accretion disk in a central region of quasar is beyond any doubts, as well as the fact that the r eal appearance of 



the qu asar core can be quite different from very simplified theoretical models in question. On the other hand. lMortonson et al 



( 20051 ) argue that the accretion disk can be modelled with any brightness profile (Gaussian, uniform, etc.), and this model will 
agree with the available data provided that an appropriate source size is chosen. Anyway, in these considerations, different 
effects come into play, and the consistent treatment requires to consider all the corrections to the amplification including those 
obtained in the present paper. 
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APPENDIX A: ANALYTICAL EXPANSION METHOD 



Here, we present a justification of the method of analytical expansions near the folds used in Section l371] following: IAlexandrov et al 
I 2003h . This problem is not trivial, because the Jacobian of the lens mapping on the fold is equal to zero. 

We say that an ana lytic function f{t) has order fc at t = , if /''^O) = for i = 0, l,...,fe - 1, and /''''(O) / 
I Poston fc Stewartlll97i '): then we write f(t) = 0{t''). 

It is well known that the mapping 



y = F{x) 



(Al) 



of two-dimensional manifolds (e.g., mapping ([3])) in a neighborhood of the fold (x = y = Q) can be reduced by the coordinate 
transformations 



u^u{x), v = v{y), x = x{u), y = y{v) 

to the normal form (see, e.g., Poston & Stewart 19781 : Petters et al. 2001) 

U ^ V : V-i — U\ , V2 — {U2f . 



(A2) 



(A3) 



We suppose the initial mapping (|A1[) is analytic; then transformations (|A2[) are also analytic. 

Let the source move along a parameterized curve v (t) such that (i) 112 (t) > for t 0, v (0) — 0; (ii) the functions 
Vi (t) are analytic at t = 0; (iii) the order of V2 (t) is 2k for some integer k ^ 1. Conditions (ii) and (iii) mean that V2 {t) = 
t^'' {q + X (^)) 1 9 7^ 0, and x (i) is analytic function at f = so that x (0) = 0. In view of (i), we have q + xi't) > 0. 

We have two obvious solutions of Eqs. (|A3|) : 



U(±)(f) = {Vi{t),±t''y/q + x{t)}, 



(A4) 
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which, in view of g > 0, are analytic as functions of t at the common point t = 0. We can say that the trajectory of the source 
consists of two branches, corresponding to values t < and t > 0. Each branch has two images lying on different sides of 
the critical curve. When k is even (odd), the images which lie on the one side (on different sides) continue analytically one 
another. Thus, a suitable choice of the source trajectory v{t) leads to the analyticity of its critical images as functions of the 
parameter. 

We now intend to reformulate the sufficient conditions of the existence of analytical solutions in the original coordinates. 
The coordinate transformation y — > u in the source plane can be written as 

vi = Ayi + By2 + -l-i (yi, ^2) , (A5) 
f2 ^Cyi+Dya + $2 (1/1,1/2), (A6) 
AD-BC^ 0. 

Here, $i (3/1,2/2) are analytic functions having the Taylor expansions starting from the second order. 

Then we assume that the terms linear in Xi are present only in the first component, F\{x), of the lens mapping (IA1[I . 
Then, in Eq. (IA6p . C = 0. This statement follows from the substitution of in (|A6|l . the application of the transformation 
X = X (u), and the comparison with the second equation in (IA3p . 

Suppose that 

(a) the curve y(t) lies on the inner side of the caustic (except y{Q) — 0), 

(b) y{t) is analytic and such that yi{t) = 0{t"^),y2{t) = ©(t^*) ; m > fc ^ 1, and m and k are integers. 
Then Eq. IfAffil yields V2{t) = 0{t'^'') and conditions (i),(ii), and (Hi) are fulfilled. 

In summary, we obtain the following statement. 

Let the function F{x) of the lens equation (Ailj be an analytic function of x in a neighborhood of the fold critical point 
X — 0, let y — be the corresponding caustic point, and let F2{x) do not contain linear terms. Let also y{t) be an analytic 
vector-function which represents the curve satisfying conditions (a),(b). 

Then the lens eauation \Al\ has two solutions X(^j^){t) and a;(_)(i) which are analytic in t and represent the critical images 
of the curve y{t) . 

These conditions are sufficient ones, and they do not exhaust possible combinations of k, m. Typically, the statement 
is true also when m — k. This statement allows us to use analytical expansions along test source trajectories to obtain the 
critical solutions near the folds. In the main body of this paper, we have used the family of straight-line trajectories of the 
source Viit) = ajt^ (k = 1, m = 2) in a neighborhood of the caustic on one side from the tangent to the caustic. In work 
I Alexandrov et allboosh . we also considered the version with parabolic curves yi (t) = at, 2/2 (t) = bt^, which allowed us to 



find the solutions in the approximation of a parabolic caustic. 



APPENDIX B: ANALYTICAL ITERATION METHOD 

In this Appendix, we show that the critical solutions of the lens equation (|3]) (after the substitution of Q) can be represented 
in the form the functions p,r,s,w being polynomials in t,yi, and 2/2 on every step of the approximation procedure. 

Ifere, we use the notations which are independent of the rest of the paper. After a simple change of the variables 
{yi ~ {by I — ay2)/2b, 1/2 — —y2/b), Eq. ([6j can be written in the form 

yi = xi + t''^a„,m{t)xix^ , 

ri,m 

y2 = xl + t^'^b„^,n(t)xiX^ , (Bl) 

where the indices m,n take on integer nonnegative values; t is a small parameter; the coefficients a„^rn{t) and bn,,n{t) are 
finite-order polynomials in t. In fact, all considerations below can be performed in the case where the r.h.s. of (|B1[I is an 
analytic functions in a;i,a;2, and the small parameter t. However, in this paper, we deal with finite-order approximations. 
We search for a solution in the form 

xi=p + t£r^^, X2 = ts + e\i^; e = ±1. (B2) 

After the substitution of Eq. (|B2} in HB1|I . we separate the terms containing integer and half-integer powers of w, e.g., 

a„^m{p + £tr^)"{ts + sy/w)"^ = Ao + ew^^^Ai, 

where 
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Ao=J2an,m J2 C'^C!C{tr)''p"-''{tsr-'''w'^, if = is an integer, 

k,k' 
k + k' = 2K 

Ai=^a„,„ J2 C'^C^{tr)''p"-\tsr-'''w'^, R = is an integer. (B3) 

k,k' 
k + k'=2K + l 

Here, are the binomial coefficients, = for k > n and for fc < 0; and k, k' are nonnegative integers. 

Then we equate separately the terms with integer and half-integer powers of w on both sides of the relation following 
from the first equation of system (|B1|> and obtain 

p = yi+tP{t,p,r,s,w), (B4) 
where 

P{t,p,r,s,w) = -Ao, 
and 

r + Ai=0. (B5) 

Analogously, the second equation of (|B1|) yields 
w = y2 + tW{t,p,r, 3,10), (B6) 
where 

W{t,p,r,s,w) = ~ts^ — Bo, 
and 

2s + Bi = 0. (B7) 



Here, Bo and Bi are defined similarly to Ao and Ai by Eq. (|B3|) with the replacement of the coefficients, a„,m — >■ 6„,m. 

Separating the terms of the zero order with respect to r and s (or with respect to t) in Eq. (IB3p . we rewrite Eq. (|B5|l as 
follows: 

r + ^ an,2K+ip"w'^ +tU{t,p,r,s,w) = 0. (B8) 
Here, 

U{t,p,r,s,'w) = 2_^an^m C„Crnt rp s w . 

1,™ fc.fc' 

k + k'=2K+l 
k + m-k' >0 



Substituting p and w from Eqs. (IB4I) and (|B6P in Eq. (|B8|I . we get 

+ ^an,2K+i(yi +tP)"(y2 + tt/(t,P,r,s,^j;) = 0. 



r 



This relation can be written as 

r ^ f{t,yi,y2) + tR{t,p,r,s,w), (B9) 
where 

/(i.yi.yz) = - '^an,2K+iyiV2 , 

R{t,p,r,s,w) = -^^a„,2K+i ^ C^Cxj/" "''yf "'^ P*" t*"''''' - [7(i,p, r, s, w). 

fc,fc' 
fc + fc'>0 



An analogous consideration of Eq. (|B7P yields 
s + i ^ bn,2K+ip"'w'^ + ^V(t,p,r, s,w) = 0, 



2 



where 
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n,m k,k' 

k + k'=2K+l 
fe + m — fc'>0 

Then we substitute p and w from Eqs. (|B4|l and (|B6|l to obtain 

s = g{t, yi, j/2) + ^^(t,^, r, s, w), 
where 



(BIO) 



n X 



S{t,p,r,s,w) 



fc + fc'>Q 



Thus, we have the system of equations (|B4|I . (|B9p . (|B10p . (|B6|) for p, r, s, and w that can be represented in the form 



X = F{t,Y)+tG{t,X,Y), 
where 



(Bll) 



P 
r 
s 
w 



yi 
y2 



Fit,Y) 



yi 

fit,yi,y2) 
9{t,yi,y2) 
y2 



Git,X,Y) 



P{t,p,r,s,w) 
R{t,p,r,s,w) 
S{t,p,r,s,w) 
W(t,p,r,s,w) 



System (IBllf) is ready for iterations 
X(„) =F(f,y)+tG(t,X(„„i),V), n = l,2,..., X^o)=F{t,Y). 



For a sufficiently small t, the iteration process converges due to the contraction mapping theorem. It is worth to note 
that, at every iteration step, we obtain an approximate solution in the form of finite-order polynomials in Y. This is obvious 
from the explicit form of the functions f{t, yi, 1)2), g(t, yi, 7/2), P{t,P, r, s, 10), R{t,p, r, s, w), S{t,p, r, s, w), and W{t,p, r, s, w). 
The application of the above procedure to system ^ yields solution pH14l) . 



APPENDIX C: AMPLIFICATION FOR A POWER-LAW SOURCE 



Here, we use Eq. (|18|) to derive the magnification of an extended centrally symmetric source with a power-law brightness 
distribution. Two different types of "power-law" distributions can be found in the literature on the gravitational lensing. In 
particular, there are the distributions (|Shalvapinii2001 : Shalvapin et al..,200^ ') 

where p > 1 is the power index, the source centre is at the coordinate origin, the distribution HC1|) is normalized to unity, and 
L is related to the r.m.s. radius Rrms as = {p — 2) Rrmsy 

Rrms = J dyidy2y^Ip~\y). 

For fixed Rrms and p — > 00, the brightness distribution (|C1I) tends to the Gaussian one. 

Along with (|C1|) . the models for limb darkening are also often used in microlensing studies (see, e.g.. lDominildl2004l ) 

4'''(y) = ^S(|yl/L;g), E{^- q) = 0(1 ^ ^ fy , (C2) 



where L stands for the source radius, and L'^ = {q + 2) Rrms- Here, we assume q > 0. Linear combinations of distributions 
(|20|l . HCip . and HC2|) with different parameters yield rather a wide class of symmetric source models. 

For brightness profile (|C1[) . the total microlensed fiux (|18p that describes a variable contribution of the critical images, as 
the source crosses a fold caustic, is the convolution of (|C1[) with (1191) . The result for the amplification factor involves integrals 
that can be expressed via the hypergeometric function 2F1 (|Bateman fc Erdelvi.,1953i l: 
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r(p- 



for fc = 0, 1, -B (a::,2/) being the Beta-function. 

We extend ()C3|) to fc = —1 having in mind the definition of (y)^^^^ , so that 



(C3) 



(C4) 



Like in Subsection 14.31 we introduce the normafized coordinates of the source center s = Yi/L, h = Y2/L. Now, the 
amphfication due to critical images takes on the form 



Kp {s,h) 



2y/-K\b\L 



*0,p(^) +^ 



(C5) 



The zeroth approximation to this formula has been derived bv lShalvapirJ (|200ll ') . 

In the case of the model with limb darkening (|C2p . the critical images disappear when the source lies on the outer side 
of the caustic (i.e., for h < —1). The substitution of (IC2I) and p9l) in p8[) yields the total amplification of critical images as 



2y/TT\b\L 

where we denote 



Xo,, (h) + L 



" X_i,,+i {h) + QsXo,, [h) - ^ s'X_i,, [h) 



8 (g + 2) 



(C6) 



Xfe., (ft) 



r(g + 2) 



/-^ H(y-/i;q + l/2)dt/, fc = l,2. 



u 

Like the previous analogous cases, we define 
X_i,, {h) = 4 (g + 1) (ftXo,,-i - Xi,,_i) . 
We have 



Xfc,, {h) = 2"+^ (1 + hY+''+ 



,r (g + 2)r(fc + i) 

r(g + fc + 2) ■ 



Fi -g- -,g+ -;g + A; + 2: 



l + h 



for — 1 < /i < 1 and 
Xfc,,(ft) = ^(ft + l)* 



Fi(g+|,i-fc;2g + 3; 



h+1 



for h> I. 

The functions '^o,p, '^i.p, ^'-i.p, and Xo,q, Xi,^, X_i,q are analogous to $0, 4$2 from (25), respectively; they have a 
similar qualitative behavior and the same asymptotics. Note that all formulas for the extended source models are transformed 
into that for point-source amplification l|19p . as the source size tends to zero. 



APPENDIX D: EXPLICIT COMPARISON OF THE FIRST-ORDER FORMULAS WITH 
COUNTERPARTS FROM THE PAPER BY KEETON, GAUDI & PETTERS (2005) 

Here we compare our expressions for the first order corrections with that of Appendix A2 from ( Keeton et al. 20051 ). further 
KGP. This is especially relevant because KGP considers the general case of the lens equation without supposition on harmonic 
potential. 

Initial equations (A7-A8) of the lens mapping of KGP are as follows 

= KOi - (366)1 + 2/e»ie»2 + gel) - {Akel + smejez + 2nei0l + pel) , 
= - {.fej + 2ge»i6'2 + shel) - (mel + 2nele2 + Speiel + 4rel) . 

These equations need to be compared to our Eqs. (|3]|6]). The correspondence between the coordinate notations is -O Xi, Ui -f-)- 

In the general case the lens potential <I> (x) obeys the equation 

A$ = 2cr(a;), (Dl) 

where a (x) is the normalized surface density of the continuous matter. Coordinates are chosen to diagonalize the symmetric 
matrix Aij — dyj/dxj] r^- Its eigenvalues ar e Ai = 1 — cr(0) -I- r(0) and A2 = 1 — cr(0) — r(0), where F is modulus of the 
complex shear (see, e.g.. lSchneider et al.lll992l ). At the critical point one of the eigenvalues is zero (the coordinates are chosen 



so that dy2/dx2\ 



0), in this case we have dyi/dx\\ 



_ft' = 2[l-a(0)]. 
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Now we consider the correspondence of the expansion coefficients. For example Eq. ([T]) yields 3e = |$,iii; g — |$,i22- 
Under the assumption that during HAE a = ao — const we have 3e — —g — >■ —a. Analogously / = — 3/i — )■ — fe, 
4fc = — |n = 4r — > —c,p = —m — >■ d. 

For critical solutions near a fold caustic, KGP found Eqs. (Af4-A15) which we repeat below along with our analogs 



± Shui - gu2 , ,3/2 ^ - ,2 byx ~ ay2 2 , ^ /,3\ 



'—^ + O{0^'^ ^ (D2) 



L-y2 abyi- [a -c(l^ao) y2 
X2=X2t = ^e^—t 2bHl^ao) ' + ^ ) ' 

For ao = these expressions are equivalent to xi = iio + O (t), £2 ~ £20 + X2\t + O (t^), where iio, 520, i2i are given by 
(|7I8|) . Thus (from the viewpoint of our approach) KGP derives the first coordinate in zero approximation. However, as KGP 
point out this is sufficient to find the first order correction for the amplification. Namely, for the corresponding Jacobian of 
the lens mapping, KGP found (Af7): 

(m*) = ±2K^/^3hu2^^^^ + 4 ~ ^^'^ ^ ^^'^^ "'^ ^ ^ ^ 

J = 4et (1 - ao) y^^b^ + ^t^ [a^ + b^-c{l- ao)] §2 + (t^) . (D3) 
For ao = the latter equation is the same (in the first approximation) as Ea. (|f5p . 

For ao 7^ the approximate solutions follow from the corresponding ones with ao = by means of the rescaling 
-5> ?/i/(f - ao), a-^a/(l-ao), 6->6/(f-ao), c ^ c/(f - ao), ... (D4) 
Analogous rescaling for the Jacobian is J — > J (f — ao)^; it is this rescaling of the variables has been mentioned in Section[2l 
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